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Abstract. —A Habitat Suitability Index Model (HSI), following the U.S. Fish and Wildlife Service HSI Model 
Series, is described for the black-footed ferret. The literature on which the model is based is reviewed, and model 
assumptions and structure are discussed. A realistic model is specified with variables and their functions that embody 
the critical spatial and resource heterogeneity characteristic of the broad geographic environment ferrets occupy. It 
assumes that ferrets can meet year-round habitat requirements within prairie dog colonies providing: (I) prairie dog 
colonies are large enough, (2) burrows are numerous enough, and (3) adequate numbers of prairie dogs and alternate 
prey are available. Five habitat variables are identified: VI is the frequency distribution of colony sizes, V2 is the total 
area of colonies, V3 is burrow opening density, V4 is intercolony distance, and V5 is prairie dog density. Variables are 
compensatory. As more data become available and our understanding of ferrets expands, the basic model design can 
readily incorporate improvements without radical restructuring. 


Habitat models are an attempt to describe and 
quantify an animal’s essential habitat require¬ 
ments or “life requisites” and are therefore a 
useful tool in habitat evaluation. The Habitat 
Suitability Index (HSI) Model Series, developed 
by the U.S. Fish and Wildlife Service (USFWS), 
provides habitat descriptions for several species. 
These models are useful for assessment of im¬ 
pacts on wildlife and habitat management (US¬ 
FWS 1980a, b) and may prove especially valu¬ 
able in endangered species management, where 
determination of habitat quality and suitability is 
often critical for management and continuation 
of the species. HSI “models should be viewed as 
hypotheses of species-habitat relationships ra¬ 
ther than statements of proven cause and effect 
relationships” (Schamberger et al. 1982:1). 

This paper applies the HSI Model format to 
the Meeteetse, Wyoming environment of the 
black-footed ferret (Mustela nigripes; BFF) as 
generally described by Clark et al. (Descrip¬ 
tion and history, 1986) and more specifically 
by Forrest et al. (1985)(Fig. 1). Applications 
and uses of the model are: (1) to compare other 
areas to BFF habitat at Meeteetse, (2) to use 
those comparisons to select areas to be 
searched for BFFs, and (3) to select suitable 
areas for transplant sites. Our use of the HSI 
format closely follows the USFWS (1981) and 


parallels applications by Allen (1982a, b, 
1983, 1984) for other species. 

Our use of the HSI model for BFFs incorpo¬ 
rates several recent improvements on the 
roles of ecological models: (1) We stress model 
reality of a single species more than focus 
upon model precision or generality (see 
Levins 1966, Rosen 1978, Kaiser 1979, Pielou 
1981). (2) Few highly measurable variables 
dictate the HSI, and, although some are colin- 
ear, together they contain high explanatory 
power, at the same time allowing comprehen¬ 
sible results and simplified sensitivity analy¬ 
sis. This reflects the growing consensus that 
there is no apparent relation between model 
complexity and predictive utility in any field 
of forecasting (e.g., Ascher 1978, K. E. F. 
Watt personal communication). (3) Our model 
uses nonlinear representations of variables, 
rather than linear, because those more accu¬ 
rately express the dynamic nature of biologi¬ 
cal responses and realistic species-habitat re¬ 
lations (Whittaker 1975, Green 1979, West- 
man 1980, Johnson 1981, Meents et al. 1983). 
Nonlinearity permits us to mimic more realis¬ 
tic biological processes that involve thresh¬ 
olds and limits and the smoothed transitions 
between them (Holling 1985, J. R. Krebs per¬ 
sonal communication). (4) The model vari¬ 
ables and their functions embody the critical 
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Fig. 1. Photographs of black-footed ferret habitat (prairie dog colonies and prairie dogs) and ferret predation. Photos 
by Tim Clark. 

A. White-tailed prairie dog colony occupied by ferrets. 



B. Black-footed ferret at prairie dog burrow. 
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D. Black-footed ferret with prairie dog prey. 
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importance of spatial and resource heterogene¬ 
ity. The structural simplicity of the BFF-prairie 
dog (Cynomys spp.) community promotes a de¬ 
sign where all variables directly assess spatial 
patchiness and resource variability, consider¬ 
ations that have pivotal impact on population 
dynamics and population viability (reviews in 
Steele 1974, Wiens 1974, Southwood 1977, 
Shugart 1981). 

The outeome of the above four features is only 
a slight increase in model complexity traded for a 
dramatic increase in ecological reality. Perhaps 
of equal benefit is the ease of model validation. 
As more data become available and our under¬ 
standing of BFFs expands, the basic model de¬ 
sign can readily incorporate improvements with¬ 
out radical restructuring. Data sets already 
completed and cited below could likely be 
reevaluated with future model versions. 

This HSI application for the BFF draws on 
Clark et al. ( Description and history , 1986) and 
Forrest et al. (1985), who describe the Mee- 
teetse, Wyoming, BFF study area (1981-1985) 
and its use by BFFs as well as all the data from 
the Mellette County, South Dakota, BFF study 
(1964-1974). Because of the localized nature and 
limited size of these two study areas, this HSI 
model will likely require updating if BFFs are 
found in other areas in different ecological set¬ 
tings. In the meantime, this HSI model can 
serve as a useful tool in BFF recovery planning 
to evaluate proposed transplant/relocation sites. 

Background 

Requests for evaluation of BFF habitat have 
been frequently mentioned in the literature. 
The Black-footed Ferret Recovery Team 
(1978) requested research to define compo¬ 
nents of a prairie dog colony necessary to sup¬ 
port BF'Fs. The BFF Recovery Plan al so notes 
the need to establish ideal habitat sites for 
successful introduction of transplanted BFFs 
(see Linder et al. 1972) . The South Dakota 
BFF and Prairie Dog Workshop in 1973 sug¬ 
gested several BFF management needs, in¬ 
cluding a definition of habitat (Hillman and 
Linder 1973, Stuart 1973, Erickson 1973). 
Others have discussed the need for BFF pre¬ 
serves and habitat descriptions (Clark 1976, 
1984, 1986). Flath and Clark (1986) described 
historic prairie dog distributions in Montana 


for the period 1908-1914. This early Montana 
situation probably represented a habitat setting 
in which BFFs evolved among the complex in¬ 
terrelationships of species and environmental 
interactions of the prairie dog ecosystem. 

Hillman et al. (1979) described prairie dog dis¬ 
tribution in the area occupied by BFFs in South 
Dakota. Their description was widely used by 
management agencies as a guide to the number 
and spacing of prairie dog colonies to be left after 
prairie dog eradication programs. 

Clark et al. ( Description and history, 1986) 
provided a descriptive and historical overview of 
the Meeteetse BFF environment. Forrest et al. 
(1985) noted that BFFs are restricted to a prairie 
dog complex—a group of prairie dog colonies 
distributed so that individual BFFs can migrate 
among them commonly and frequently. The 37 
colonies of the xMeeteetse complex (total size 
2995 ha) were described and their occupation 
history by BFFs noted. The average density of 
adult BFFs was 1 BFF/56.6 ha. Burrow open¬ 
ings, based on literature reviews, are correlated 
with the number of prairie dogs present 
(r 0.71). High burrow densities are desirable 
for BFFs in that they provide added protection 
from predators and shelter from the elements. 
Colonies greater than 100 ha supported more 
than two resident adult BFFs, whereas colonies 
from 12.5 ha to 102.0 ha supported only one BFF 
throughout the year. BFFs traveled among the 
colonies, but to an unknown extent. BFFs may 
use burrows at low densities and colonies of 
small size in travels between larger colonies. 
BFFs moving between colonies have a greater 
chance of finding another colony if the colonies 
are large and close together. 

Several bibliographies of BFFs (Harvey 1970, 
Snow 1972, Hillman and Clark 1980, Casey et al. 
1986) and of prairie dogs (Clark 1971, in prepara¬ 
tion, Hassien 1973) exist. These also serve as 
background for this HSI model. General infor¬ 
mation on BFFs is summarized in the bibliogra¬ 
phies listed above, in primary sources from 
South Dakota studies (e.g., Hillman 1968, Hen¬ 
derson et al. 1969, Fortenbery 1972), and, more 
recently, from Meeteetse, Wyoming (e.g., 
Clark et al., Description and history, 1986; 
Clark et al.. Descriptive ethology, 1986; Camp¬ 
bell et al. 1985, Richardson et al. 1985; Forrest 
et al. 1985, Biggins et al. 1985). 
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Habitat Use Information 
Overview 

A member of the family Mustelidae, the 
BFF is the only ferret native to North America 
(Hall 1981) and is perhaps the rarest and most 
endangered mammal species on this conti¬ 
nent (Cahalane 1954, Hillman and Clark 
1980). BFFs are solitary except during breed¬ 
ing and maternal care of young and are pri¬ 
marily nocturnal. They prey on prairie dogs, 
whose burrows they also use for cover and 
litter rearing. 

Food 

The BFF relies on prairie dogs as its primary 
food source, although other prey, both live 
and dead, are taken in considerably lesser 
amounts (Hillman 1968, Henderson et al. 
1969, Sheets and Linder 1969, Sheets et al. 
1972, Clark et al. 1985). Sheets et al. (1972) 
found 91% of 82 BFF scats from South Dakota 
contained prairie dog remains, and Campbell 
et al. (unpublished data) found 87% of 86 BFF 
scats from Meeteetse contained prairie dog 
remains. Prairie dogs, on this basis, compose 
the major BFF food. 

Stromberg et al. (1983) generated a preda¬ 
tor-prey model of metabolizable energy re¬ 
quirements that estimated: (1) annual prey 
requirements for one reproductive female 
BFF and her litter of four and (2) prairie dog 
population sizes needed per BFF. Powell et 
al. (in press) estimated BFF winter energy 
expenditure (about 104 kcal/day) and prey re¬ 
quirements (about 20 prairie dogs from De¬ 
cember through March) at Meeteetse. A lac- 
tating female with four young are predicted to 
need six times the winter estimate, or about 
one prairie dog per day in summer. 

Water 

BFFs apparently satisfy water requirements 
through prey consumption and have never 
been observed in the wild drinking free wa¬ 
ter. Henderson et al. (1969) reported that 
captive BFFs drank water irregularly. L. 
Richardson (unpublished data) watched a 
BFF eating snow at Meeteetse. 

Cover 

Cover for BFFs is provided by prairie dog 
burrows, which are used for predator avoid¬ 


ance and thermal cover throughout the year 
(Clark et al. 1985, Richardson et al. in press) 
Any prairie dog burrow is assumed to be suffi¬ 
cient to satisfy BFF cover requirements. 
Higher burrow densities provide greater 
cover. 

Reproduction 

Reproductive habitat requirements for 
BFFs are assumed to be identical to food and 
cover requirements described above because 
all BFF activities are associated with prairie 
dog burrow systems (Clark et al. Descriptive 
ethology , 1986; Richardson et al. in press; 
Forrest et al. 1985) . Large, mounded, multi- 
entranced burrows may be important for litter 
rearing because of their presumed extensive 
tunnel network. 

Interspersion 

A picture of BFF home range patterns is 
emerging from research efforts at Meeteetse. 
A single adult male’s range may encompass 
home ranges of several females, which show 
much smaller ranges (Richardson et al. un¬ 
published data). Females remain with their 
litters until late summer, when young become 
independent (Henderson et al. 1969, Clark et 
al. Descriptive ethology, 1986). BFFs appear 
to have a typical mustelid spacing pattern de¬ 
scribed by Powell (1979), Forrest et al. (1985), 
and Richardson et al. (in press). More infor¬ 
mation is needed on BFF home ranges and 
movements, dispersal of young or adults, and 
inter- and intrasexual interactions. 

Interspersion characteristics of BFFs repre¬ 
sent a two-dimensional management consid¬ 
eration—individual and populational. Indi¬ 
vidual interspersion patterns are better 
known than populational interspersion pat¬ 
terns required for minimum population sizes. 
A resident female snow-tracked from Decem¬ 
ber through March used 16.0 ha and was over¬ 
lapped by a resident male that used 136.6 ha 
(Forrest et al. 1985). Studies of radio-collared 
BFFs show a young female used 12.6 ha in 
October and November (Biggins et al. 1985). 
Population interspersion is dependent on the 
size, configuration, and intercolony distance 
of prairie dog colonies making up the com¬ 
plex. Data show that, if colonies are too small 
and intercolony distances are too large, then 
BFF populations cannot sustain themselves. 
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The search for food (energetics) becomes pro¬ 
hibitive, avoidance of predators becomes dif¬ 
ficult or impossible, and adequate thermal 
cover is rare or nonexistent, all reducing both 
individual and population survival. 

Special Considerations 

Successful management of BFFs depends 
on maintaining adequate numbers and areas 
of prairie dog colonies. Minimum viable pop¬ 
ulation (MVP) sizes and area requirements for 
BFFs were addressed by Groves and Clark 
(1986). Additional estimates of these variables 
are underway by Shaffer et al. (in prepara¬ 
tion), who are modeling effects of both demo¬ 
graphic and environmental stochasticty on 
BFF populations of varying sizes. The MVP 
represents a threshold below which popula¬ 
tions are not self-sustaining. Populations may 
persist for a long time below the MVP, but 
probably at a loss of adaptability and a high 
susceptibility to local extinction. Groves and 
Clark (1986) noted that the genetic method of 
determining MVP for the Meeteetse BFFs 
estimated that about 200 animals are needed 
for maintenance of short-term fitness. The es¬ 
timated 200 animals needed is about four 
times the number of breeding adults esti¬ 
mated to currently exist at Meeteetse (Clark 
1986). 

Poisoning and shooting of prairie dogs 
should be prohibited from areas where BFFs 
occur as well as from other selected portions of 
prairie dog range. Hubbard and Schmitt 
(1984) suggested a “refugia” concept of man¬ 
aging prairie dogs in which relatively large 
areas are omitted from poisoning and other 
disturbance. They suggested that refugia be 
large enough to support a BFF MVP and 
based such area estimates on the Stromberg et 
al. (1983) predator-prey model. Clark (1986) 
outlined a series of management guidelines 
for BFFs. 

Differences in black-tailed (C. ludovi- 
cianus) and white-tailed prairie dog colonies 
have been noted (Tileston and Lechleitner 
1966, Campbell and Clark 1981, Clark et al. 
1982). Black-tailed colonies often show great¬ 
er prairie dog and burrow opening densities— 
two important variables of BFF habitat. Satis¬ 
fying habitat requirements for BFFs on 
white-tailed colonies as described in our HS1 
model is assumed also to satisfy habitat re¬ 


quirements on black-tailed and Gunnison’s 
(C. gunnisoni) prairie dog colonies. 

Application of Habitat Suitability Model 

Model Applicability 

Geographic area .—Although this model 
was developed on data from the only two BFF 
populations ever studied, it should apply 
throughout the historic range of the BFF until 
additional BFF populations in different eco¬ 
logical settings are found, studied, and results 
show it does not apply. Even though a single 
prairie dog colony cannot support a BFF MVP 
(unless it is extremely large), it can potentially 
support one or more individuals. Therefore, 
any prairie dog colony should be considered 
potential BFF habitat. Historic and current 
land use patterns affect the quality of BFF 
habitat. A constellation of prairie dog 
colonies, described by Clark et al. ( Descrip¬ 
tion and history, 1986) and Forrest et al. 
(1985) as a prairie dog “complex, ’’ is needed to 
support a BFF MVP. 

Season .—This model has been developed 
to compare year-round BFF habitat at Mee¬ 
teetse to habitat in other areas. Because 
prairie dogs may become torpid or hibernate 
over winter at northern latitudes, it is recom¬ 
mended that evaluation take place when 
prairie dogs are active and when snow cover is 
minimal or absent: late May to late June is 
recommended. 

Cover Types .—This model compares the 
BFF habitat at Meeteetse to other potential 
BFF habitat in all cover types where prairie 
dogs are found. 

Minimum Habitat Area .—Minimum habi¬ 
tat area, as discussed for BFFs by Forrest et 
al. (1985), is defined as the amount of contigu¬ 
ous habitat that is required before an area will 
be occupied by a species (Allen 1982a). We 
recommend that a preliminary estimate of 
4,000-6,000 ha of prairie dogs is needed to 
support a MVP of 100 BFFs (Forrest et al. 
1985, Groves and Clark 1986). 

Model Review. —Drafts of this model were 
reviewed by our colleagues in the Idaho State 
University/Biota Research and Consulting, 
Inc. ferret study team-Steven Forrest, Louise 
Richardson, Tom Campbell, and Denise 
Casey; Arthur Allen, Habitat Evaluation Pro¬ 
cedures Croup, USFWS; Wayne Brewester 
and Ronald Crete, Office of Endangered Spe- 
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Fig. 2. The relationship of habitat variables, life requisites, and cover types to the HS1 for the black-footed ferret. 


cies, USFWS; Donald Streubel, Department 
of Biological Sciences, Idaho State University; 
Craig Groves, Idaho Heritage Program, The 
Nature Conservancy; Mark Stromberg, The 
National Audubon Society; John Hubbard, 
Endangered Species Program, New Mexico 
Game and Fish; John Cadaand Dennis Flath, 
Nongame Program, Montana Department of 
Fish, Wildlife, and Parks; Harry Harju, Wyo¬ 
ming Game and Fish Department; and Sid 
England and Dale Lott, Department of 
Wildlife and Fisheries Biology, University of 
California-Davis. Improvements and modifi¬ 
cations suggested by these persons are appre¬ 
ciated and were incorporated into this model. 

Model Description 

Overview .—The BFF can meet its year- 
round habitat requirements within prairie 
dog colonies providing: (1) prairie dog 
colonies are large enough, (2) burrows are 
numerous enough, and (3) adequate numbers 
of prairie dogs and alternate prey are avail¬ 
able. This model therefore assumes that re¬ 
producing populations of BFFs use only 
prairie dog colonies, and habitat evaluation 
based on this model considers only the life 
requisites provided by such colonies. BFFs 
do not rely solely on prairie dogs for food, but 
breeding populations may depend on prairie 
dog colonies with their host of associated ver¬ 
tebrates, many of which are known food 
items. It assumes that these colonies will 
provide a sufficient prey base (including alter¬ 
native prey) and sufficient burrow openings 
for predator evasion and as sites of litter rear¬ 
ing, thus providing maximum potential for 
BFF habitat. Ecological differences in habitat 


may be found if future populations of BFFs 
are discovered, or if BFFs are found on areas 
other than prairie dog colonies. 

The following section documents the logic 
and assumptions used to translate habitat in¬ 
formation for the BFF to the variables and 
equations used in the HSI model. Specifi¬ 
cally, this section covers: (1) identification of 
variables used in the model, (2) definition and 
justification of the suitability levels of each 
variable, and (3) description of the assumed 
relationship between variables. The BFF 
habitat variables have been grouped into two 
sets: (1) an aggregrated set of four variables 
that assess cover/reproduction as life requi¬ 
sites and (2) a single life requisite variable for 
food. Figure 2 illustrates the relationship of 
habitat variables, life requisites, and cover 
type for the BFF. The five habitat variables 
identified under the two life requisite cate¬ 
gories are: VI is the frequency distribution of 
colony sizes, V2 is the total area of colonies, 
V3 is burrow opening density: average num¬ 
ber of burrow openings per ha of colony, V4 is 
intercolony distance: mean distance between 
colonies (these four variables are grouped un¬ 
der the cover/reproduction life requisite), and 
V5 is prairie dog density: mean number of 
prairie dogs per ha (this variable is the food life 
requisite). The aggregrated variables are 
viewed as compensatory (i.e., an increase in 
one variable will increase the HSI, but not the 
suitability of other variables) and thus are 
combined to produce a single HSI. The limit¬ 
ing factor method is suggested for evaluating 
resulting values of the two variable sets. 

Cover!reproductive component. —BFFs re¬ 
ly on prairie dog burrows for cover and litter 
rearing. Four variables are defined. 
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Variable 1 examines the relationship be¬ 
tween the distribution of prairie dog colony 
sizes in a region and its Suitability Index. 
Prairie dog colonies present at the turn of the 
century represented extremely large areas of 
contiguous prairie dog distribution (e.g., in 
Montana see Flath and Clark 1986). Such ar¬ 
eas represented a 100% prairie dog occupancy 
and were assumed to be optimal habitat for 
BFFs. By comparison more recently, Mel¬ 
lette County, South Dakota, showed 1.7% of 
its area occupied by prairie dogs, with a mean 
colony size of about 9 ha (Hillman et al. 1979). 
The Big Horn Basin of Wyoming containing 
the Meeteetse BFFs has about 1.7% of its area 
occupied by prairie dogs in many small, low- 
density colonies, although a few exceed 1,000 
ha (Clark et al. Description and history, 
1986). Clark et al. (1982) described several 
sample areas in New Mexico that showed 
about 1% in prairie dogs, with colony sizes 
averaging 33 ha (range 10-61 ha); in Utah 
about 1.9%, with colony sizes averaging 33 ha 
(range 2-73 ha); in Wyoming on Thunder 
Basin National Grassland about 1.3% in 
prairie dogs showing a wide range in colony 
sizes; in southern Wyoming about 3.2% in 
prairie dogs, with colony sizes ranging up to 
2,500 ha; and in another area in Utah, colonies 
averaged 125 ha (range 0.2-958 ha). The total 
sizes of these areas varied, and this fact clearly 
influenced the distribution of prairie dog 
colony size located. If a line is drawn around 
the prairie dog complex at Meeteetse (least 
polygon enclosing all 50+ ha colonies) and the 
area occupied by prairie dogs inside this 
polygon is calculated (about 130 sq km), then 
about 22% of the area is occupied by prairie 
dogs. The 50 ha figure does not mean that 
smaller colonies are not important to BFFs; 
indeed the smaller colonies are used at Mee¬ 
teetse (Forrest et al. 1985). Colony size distri¬ 
bution within this area is listed in the Ap¬ 
pendix (Table 3). 

VI is a multidimensional probability estimate 
and is not graphable as are the remaining vari¬ 
ables. The Appendix describes computation of 
VI. 

Variable 2 is the total area of prairie dog 
colonies. Assuming a BFF MVP consists of 100 
breeding adults (even though Groves and Clark 
[1986], using genetic methods, estimated 200), 


then 100 colonies of 50 ha each (about 5,000 ha) 
is required to support them. It is assumed that 
greater colony area means greater sites for 
cover and reproduction for BFFs. 

Variable 3 is burrow opening density: the 
average number of burrow openings per ha of 
colony. Colonies at Meeteetse are character¬ 
ized by burrow opening densities as low as 10 
openings/ha and up to 100+ openings/ha. 
This compares with other areas ranging 
21-135/ha for black-tails, 32-57/ha for Gun- 
nisons, and 2-64/ha for other white-tails 
(Clark et al. 1982) . It is assumed that the 
greater the burrow opening density, the 
greater the cover and sites for successful rear¬ 
ing of young. 

Variable 4 is the mean of intercolony 
(nearest neighbor) distances. This variable is 
essential for cover/reproductive require¬ 
ments but is also essential for expansion of 
BFF populations and dispersal. In pristine 
times, BFFs in large colonies may have dis¬ 
persed from their natal areas to new areas 
without ever leaving the single large prairie 
dog colony. Dispersal between colonies, 
where escape cover is minimal or absent, is 
thought to expose BFFs to high rates of mor¬ 
tality. Intercolony distance at Meeteetse is 
about 0.92 km (range 0.13 to 3.70) . In South 
Dakota intercolony distance averaged 2.4 km. 
Intercolony distance for a sample of 11 Gun¬ 
nison’s colonies in New Mexico, Colorado, 
and Utah was 2.4 km and for 33 white-tailed 
colonies in Utah and Colorado was 4.9 km. In 
winter at Meeteetse BFFs in intracolony 
movements often travel 2+ km per night 
hunting. Movements up to 8 km have been 
noted during the breeding season. It is as¬ 
sumed that the smaller the intercolony dis¬ 
tance, the higher the quality of BFF habitat. 

Food component .—Food is described by a 
single variable. 

Variable 5 is prairie dog density (number/ 
ha). High densities of prairie dogs provide 
increased opportunity for BFFs to success¬ 
fully meet their energy and nutrient require¬ 
ments as well as providing alternate prey asso¬ 
ciated in prairie dog colonies colonies. 
Additionally, a high density of prairie dogs 
means an increased density of burrows, which 
is related to the previous variables as well. 
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Table 1. Equations for determining year-round life requisites for the black-footed ferret (2.0 is included as a scaling 
factor). 


Life requisites Cover type Equation 

Cover/Reproduction All cover types where prairie dog (2x VI x V2xV3x V4 ) 14 

colonies occur 

Food Same as above V5 


Variable Relationships 

Suitability of BFF habitat depends entirely 
on attributes of prairie dog colonies. VI con¬ 
verts the distribution of colony sizes (relative 
to the total colony area) into a single SI mea¬ 
sure. V2 accounts for the total area of colonies 
relative to BFF requirements and is espe¬ 
cially discriminative in the range of MVP area 
size. V3 gauges the value of colonies in terms 
of cover (burrow opening density) and, al¬ 
though it generally covaries with food (V5: 
prairie dog density), any particular case may 
be critically uncorrelated. V4 (intercolony dis¬ 
tance) appraises the effect of colony dispersion 
in reference to BFF mobility and behavior. In 
summary, VI reflects colony size distribution, 
V2 the total colony area the size distribution 
represents, V3 the cover value of the colonies, 
V4 the spatial dispersion of those colonies, 
and V5 the food value of the colonies. 

Suitability Index (SI) graphs and equations 
for habitat variables . —This section contains 
suitability index graphs and equations that 
illustrate the habitat relationships described 
in the previous section (Fig. 2). 

Equations .—Life requisite values for the 
BFF can be obtained by combining the SI 
values through the use of equations (USFWS 
1981). A description and explanation of the 
assumed relationship between variables was 
included under the Model Description, and 
the specific equations in this model were cho¬ 
sen to mimic those perceived biological rela¬ 
tionships as closely as possible. The suggested 
equation for obtaining year-round life requi¬ 
site values for the BFF are given in Table 1. 
The four cover/reproduction variables are 
multiplied by two (a scaling factor for VI) and 
aggregrated by using the geometric mean, 
GM. We necessarily use the GM because the 
quantities involved are measured on a ratio 
scale and the variables are not arithmetic se¬ 
quences but geometric. 


Variable 

VI Distribution of colony sizes, 

R 

/ 2 N, Q) 

P(AB|N„ i) - . 2 

N, X i J 

2 / 

With n, as the number of colonies of size i, the 
resulting probability increases nonlinearly 
with increase in colony sizes (numerator) rela¬ 
tive to the size of the complex (denomina- 
tor)(see Appendix for example calculations of 
this equation) 

VARIABLE: 

V2 Total area of colonies, 



TOTAL AREA OF COLONIES (HA) 

V3 Burrow opening density (mean number of 
burrow openings/ha of colony), 









108 


Great Basin Naturalist Memoirs 


No. 8 


V4 Intercolony distance (mean distance be¬ 
tween colonies), 



V5 Prairie dog density (mean number of 
prairie dogs/ha of colony), 



HSI determination .—The 1181 for the BFF 
will equal the lowest of the SI values obtained 
for either the Cover/Reproduction or Food 
life requisite. This recognizes limiting factors. 
The fact that V2 only scales to 10,000 ha of 
total colony area reflects the importance of an 
MVP area and does not mean that even 
greater-sized prairie dog complexes are not 
more desirable. The larger the complex the 
better. The largest complex sizes available 
should be selected for BFF translocations, 
and these should exceed the size of Meeteetse 
(Appendix). An HSI approaching 1.0 is ideal 
and not necessarily attainable; that is, a math¬ 
ematical ideal or extreme to compare 
against—in actuality, perfect habitat does not 
exist. 

Application of the Model 

Definitions of variables and suggested field 
measurement techniques are presented in 


Table 2. Vegetative cover types for each vari¬ 
able are those that contain prairie dog 
colonies. The Appendix contains HSI calcula¬ 
tions for Meeteetse and two other areas. 
These are presented as examples of model 
application, for ease of application of HSI to 
other areas, and for comparative purposes. 

I n terrelat ionsh ips 

Thi ee considerations from application of the 
HSI format to the Meeteetse BFF environ¬ 
ment as described in the Appendix must be 
addressed. First, the Meeteetse HSI of 0.590 
for the cover/reproductive variables is 
midrange in a HSI range of 0-1. Prairie dog 
complexes should be located that exceed the 
Meeteetse HSI and that can support large 
BFF populations well above the MVP. It is 
the low V2 (complex size or total colon) 7 area) 
that deflates the overall HSI. Second, if high 
HSI areas cannot be located that can support a 
MVP, then a series of smaller areas showing a 
lower HSI than Meeteetse will have to be 
utilized in a complex, complementary, and 
closely managed situation. Third, application 
of the HSI format to the prairie dog area in 
Mellette, South Dakota, may show that its 
HSI is well below estimated MVP require¬ 
ments. If so, this means that management for 
a minimum area and colony size pattern as 
suggested by Hillman et al. (1979) has been 
below the area needed to sustain a MVP of 
BFFs and that new recommendations are 
needed. 

Sources of Other Models 

No habitat models for the BFF were located 
in the literature except for descriptions of 
BFF habitat by Hillman et al. (1979), Strom- 
berg et al. (1983), and Forrest et al. (1985). 

Concluding Remarks 

If the prairie dog colony size distributions 
shown in Table 3 of the Appendix are typical of 
prairie dog complexes, then in terms of prairie 
dog complex area, the sum influence of most 
colonies will be less than the few very large 
ones. Distributions with this property of ag¬ 
gregation or clumping are called contagious 
and can often be modeled by generalized dis¬ 
crete distributions (reviews in Coleman 1964, 
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Table 2. Definitions of variables and suggested field measurement techniques. 


Variable definition 

VI Frequency distribution of colony sizes 
(all inhabited by prairie dogs) 

V2 Total area of colonies (all inhabited by 
prairie dogs) 

V3 Burrow opening density (number/ha of colony) 

V4 Intercolony distance (mean distance 
between colonies) 

V5 Prairie dog density (mean number of 
prairie dogs/ha of colony) 


Suggested technique 

Accurately map colony configurations, determine colony areas 
from maps; ground surveys are best, but some preliminary 
aerial surveys may be needed first 

Total area of all colonies in the study area based on accurate 
mapping determined for VI 

Walk colonies and count holes, or sample selected areas 
Measured from the edge of a mapped colony along the 
shortest distance to the next nearest colony 
Use minimal visual counts of prairie dogs active .5-2 hrs after 
sunrise on three consecutive mornings in mid-June; live 
capture-mark/recapture population estimate in mid-June 


Douglas 1979). If we view the colonies in a 
prairie dog complex distributed as a Poisson 
variate and assume the number of ha per 
colony has a highly nonrandom logarithmic 
distribution, then we may obtain the Poisson- 
logarithmic compound distribution (a type of 
negative binomial). However, the colony size 
distributions could not be fit to this distribu¬ 
tion even when larger colonies were ignored. 
This illustrates the extreme “contagion” of 
prairie dog aggregation and, consequently, 
the disproportionate effect of such large clus¬ 
ters on the outcome of any realistic model 
concerning BFF MVPs. 

In conclusion, from a statistical standpoint, 
the most effective and, therefore, most inten¬ 
sive and accurate data collection could be con¬ 
centrated in the large colonies (in fact, as a 
nonlinear function of size) with little loss of 
accuracy. In other words, the model is very 
robust to small colony exclusion. Indeed, for 
all variables except V4, which equally weights 
colony location and therefore dispersion ef¬ 
fect, small colonies could be ignored for data 
collection in MVP-sized complexes when con¬ 
sidering cost and time budgeting. 
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Appendix 

Calculations of HSI for Meeteetse and two 
Other Areas 

Variable values are most often arrived at by 
sampling techniques which produce confi¬ 
dence limits. In addition, stratified sampling 
or consideration of subsets of an area may 
produce a large range of possible variable val¬ 
ues for calculation. For ease of computations 


and sensitivity analysis, the equations for VI- 
V5 are given below and are easily computed 
on handheld calculators. 

Structure of VI 

Colony size has been stressed in terms of 
successful reproduction and energetics. VI 
appraises this important aspect of colonies by 
producing higher values for size distributions 
containing larger colonies and disproportion¬ 
ately lower values for a distribution (given the 
same area) containing smaller colonies. The 
following analogy may improve our under¬ 
standing of this. Assume we have two BFFs, A 
and B, and wished to distribute them on their 
own one-hectare plot in a number of Total 
Areas, each of which contains a different num¬ 
ber of colonies totaling a constant area. We 
drop the two BFFs randomly over these areas 
and note where they fall. In areas containing a 
few large colonies, BFFs A and B are noted to 
land more often on the same colony; in areas of 
many small colonies, A and B rarely share the 
same colony. Formally, if we divided a sample 
space (Total Area) into n subspaces (colonies) 
where 1 to k are colonies of size i and N 4 is the 
number of colonies of size i, then the probabil¬ 
ity of any two objects (BFFs A and B) co-oc- 
curring in the same subspace (colony) is 

£n,G) 

P(AB|N,, i) - i-2 

N, x i j 

2 / 

Since ( 2 ) = x - = x ~ x . then 
2! (X - 2)! 2 

2 N, (i 2 - i) 

P(AB|N„ i) =_1 Zl _ 

K \ 2 K 

2 N f x i ) - 2 N, x i 

In reduced form, solving for P(AB|N i5 i) is 
simply combining two summations. Summate 
Nj (i 2 —i) and store in memory M x . At the same 
time summate N; x i and store in memory M 2 , 
then calculate 

M! 

M 2 2 - M 2 . 
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Table 3. Colony sizes for Meeteetse (Area I), another 
prairie dog complex (Area II), and a hypothetical complex 
(Area 111). The frequency distribution is in order of colony 
size and grouped in three size class frequencies. 


Colony 

numbers 


2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 


Colony sizes in hectares 


Area 1 
.5 
.5 
.5 
1.0 

1.5 
2.0 

2.5 

2.5 
3.0 

3.5 
5.0 
6.0 

8.5 
9.0 

9.5 
11.0 

12.5 

20.5 
21.0 
24.0 

24.5 

26.5 

30.5 

31.5 

31.5 

49.5 


Area II Area 111 


3 

4 

5 

6 

7 

8 
10 
11 
12 

14 

15 

19 

20 
21 
22 


31 

34 

49 


> 0-25 


Structure of V2-V5 

Variables 2, 3, 4, and 5 are intrinsically non¬ 
linear and are each derived from the differen¬ 
tial equation dy/dt ay“ + by + c and simpli¬ 
fied to the logistic form of 

Y = (1 + ke“ r1 )- 1 

The logistic form is particularly suitable in 
describing these variables because it depicts 
two asymptotic limits (at 0 and 1 adjustable 
toward infinity) and contains an inflection 
point around which the most rapid rate 
changes occur. For example, intercolony dis¬ 
tance (V4) reflects the ability of BFFs to inter¬ 
cept life requisites upon leaving one colony 
for another. If a straight-line 10 km is as much 
as BFFs might move in a night, then that 
value is the inflection point around which crit¬ 
ical and therefore extreme shifts in the suit¬ 
ability index (SI) occur. Of course, BFFs eas¬ 
ily move from 0 to 5 km and SI values change 
little within that range. Likewise, once a BFF 
is well past the “point-of-no-return,” say 
15-20 km, SI value shifts are also small. An¬ 
other view is that the chance of intercepting 
another colony along a radius extending from 
its home colony is a quadratic function of dis¬ 
tance moved modified by the actual mobility 
and energetic characteristics of the BFF. 


27 

50.5 

67 

X 

V2: f(x) = (1 + 20 e 0009x )' 1 

x ha 

28 

51.5 

81 


V3: f(x) - (1 + 15 e 06x )- 1 

x burrows/ha 

29 

69.5 

87 


V4:f(x) 1 - (1 H- 70 e 

x km 

30 

97.5 

143 


V5:f(x) = (1 + 200 e 8x )' J 

x prairie dogs/ha 


31 


102.0 


151 


32 

153.0 

231 


33 

183.0 

257 


34 

196.5 

617 

200 

35 

211.0 

658 

2100 

36 

230.0 

671 

2200 

37 

1307.0 

2242 

2300 > 

Total area 

2990 

5496 

6800 

Total colonies 

37 

29 

4 


All colonies <1 ha are entered as 1 ha. Since 
colony sizes are often unique numbers or are 
entered that way as data, then N s is com¬ 
pletely eliminated from the calculations (see 
example). However, if each colony area is not 
estimated for some reason, then they can he 
grouped into intervals such as 0-5 ha, 5-10 
ha, etc., in which case the midpoint can be 
used (i.e.,2.5, 7.5). 


Examples of MSI Calculations 

Table 3 contains colony sizes for the Mee¬ 
teetse complex (Area I), an actual prairie dog 
complex elsewhere (Area II), and a hypotheti¬ 
cal area (Area III). Maps of these three areas 
follow (Fig. 3). Before computing VI with this 
data, it is important to understand that al¬ 
though Area I and Area II hav e different dis¬ 
tributions of absolute colony size, the distri¬ 
butions are quite similar in colony size 
relative to their total colony area. It is V2 that 
will account for the almost double total area of 
Area II. 

VI: First, calculate 

2 Nj (i 2 - i). 

i 2 
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Fig. 3. Maps of prairie clogs complexes used in examples of calculations of HSI. Area I Meeteetse, Wyoming, Area II 

another actual complex. Area III = hypothetical complex. 


Since each colony has a different area, N drops 
out and we add the following series for Table 3 
showing colony sizes for each area: 

Area I: (I.5 2 -I.5) + (2 2 - 2) + ...+ (230 2 -230) + 
(I307 2 -I307) 1,936,862 

Area II: (3 2 -3) + (4 2 -4) +...+ (671 2 -67I) + (2242 2 -2242) 
= 6,473,430 

Area III: (200 2 -200) + (2I00 2 -2I00) + (2200 3 -2200) 
(2300 2 -2300) = 14,573,200 


Note that we do not include areas of size 1; l 2 

-1 = 0 . 

Second, calculate 

which is really only the total colony area ma¬ 
nipulated as in the first calculation: 

K \2 K 

2 NjXi ) - SNjXi, 
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Area 1: 2990 2 - 2990 - 8,937,110 
Area II: 5496 2 - 5496 = 30,200,520 
Area III: 6800 2 - 6800 = 46,233,200 


We arrive at VI by dividing the first by the 
second calculations for each area: VI = P(AB) 
= .217 for Areal, .214 for Area II, and .315 for 
Area III. 

Notice the influence of Area I and II’s 
largest colony on the outcome of the first cal¬ 
culation. Area I: 1307 2 -1307 = 1,706,942 and 
for Area II: 2242 2 - 2242 = 5,024,322. If we 
were to split Area Is 1307 ha colony into two 
separate colonies, it would decrease the value 
of VI to .131. How far apart these colonies 
would be is accounted for by a simultaneous 
increase or decrease of V4. For instance, note 
how the small 200 ha colony in Area III is a 
stepping stone between the three larger 
colonies. Its critical position is reflected by a 
lower mean intercolony distance and there¬ 
fore a higher value for V4. 

Graphs of the variable equations and the 
above values follow. 


V2 Total area of colonies 



TOTAL AREA OF COLONIES (HA) 

V3 Burrow opening density (mean number of burrow 
openings/ ha of colony) 



NUMBER OF OPENINGS/HA 


Completing the other values we obtain: 


Variable Areal Area II Area 111 


VI .217 .214 .315 

V2 (ha) .424(2990) .876(5496) .958(6800) 


V3 (burrows/ha) .671(57) 
V4 (km) .980 (.92) 

HSI = .590 

V5 (dogs/ha) .214(5) 


HSI for: 

Areal - (2 x .217 x .424 
Area II = (2 X .214 X .876 
Area III = (2 x .315 X .958 


.385(37.3) .992(125) 
.980(1) .969(2) 

.613 .873 

.155(4.5) .987(12) 


x .671 x .980) 1/4 = .590 
X .385 x .980) 1/4 - .613 
x .992 x .969) 1/4 - .873 


Variable 


VI Distribution of colony sizes 



Using the above equation for VI, area I .217, Area II 
.214, and Area III .315. 


V4 Intercolony distance (mean distance between colonies) 


I and II- 0.980 



INTERCOLONY DISTANCE (KM) 

V5 Prairie dog density (mean number of prairie dogs/ha 
of colony) 



PRAIRIE DOGS/HA 








